Conventional Physics Explanations for the NuTeV sin 2 6w 
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*n ■ Abstract 

The NuTeV experiment has measured sin 2 6w^° n = 0.2277 ± 0.0013(stat) ±0.0009(syst), 
approximately 3 standard deviations above the standard model prediction. This discrepancy 
has motivated speculation that the NuTeV result may be affected significantly by neglected ex- 
perimental or theoretical effects. We examine the case for a number of proposed explanations. 
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1 Introduction 

a: 

Qh| The observable R~ 

CD . 

R - = -> VjJQ - ajv^N -» V,X) ^ (Ll) 



was first suggested by Paschos and Wolfenstein [ 1 1 to measure sin 2 9 W . 

R~ = --sin 2 9 w (1.2) 

under the assumptions of equal momentum carried by the u and d valence quarks in the tar- 
get and of equal momentum carried by the heavy quark and anti-quark seas. The Paschos- 
Wolfenstein numerator and denominator are independent of the sea quark momenta since a vq = 
a uq and a vq = a uq . R~ is more difficult to measure than the ratio of neutral current to charged 
current cross-sections, R v , primarily because the neutral current scattering of v and V yield 
identical observed final states which can only be distinguished through a priori knowledge of 
the initial state neutrino. Therefore, the measurement of R" requires separated neutrino and 
anti-neutrino beams. 

As a test of the electroweak predictions for neutrino nucleon scattering, NuTeV has per- 
formed a single-parameter fit using its R~ data to sin 2 9 W with all other parameters assumed to 



have their standard values, e.g., standard electroweak radiative corrections with p = 1. This fit 
determines 



sin 2 ^ n " shcll) = 0.22773 ± 0.00135(stat.) ± 0.00093(syst.) 

Ml „-(175 GeV) 2 
" 0-00022 x( L ) 

+ 0.00032 xln( ^ sg ; ). (1.3) 
150 GeV 

The small dependences in M top and -Mmggs result from radiative corrections as determined 
from code supplied by Bardin Q and from V6.34 of ZFITTER Q; however, it should be 
noted that these effects are small given existing constraints on the top and Higgs masses J4|. 
A fit to the precision electroweak data, excluding neutrino measurements, predicts a value of 
0.2227 ± 0.00037 [4], approximately 3cr from the NuTeV measurement. 

The experimental details and theoretical treatment of cross-sections in the NuTeV elec- 
troweak measurement are described in detail elsewhere [5|. In brief, NuTeV measures the 
experimental ratio of neutral current to charged current candidates in both a neutrino and anti- 
neutrino beam. A Monte Carlo simulation is used to express these experimental ratios in terms 
of fundamental electroweak parameters. This procedure implicitly corrects for details of the 
neutrino cross-sections and experimental backgrounds. For the measurement of sin 2 9w, the 
sensitivity arises in the v beam, and the measurement in the V beam is the control sample for 
systematic uncertainties, as suggested in the Paschos-Wolfenstein R~ of Eqn. (Il.ll ). For si- 
multaneous fits to two electroweak parameters, e.g., sin 2 9w an d p or left and right handed 
couplings, this control of systematics by the Paschos-Wolfenstein R~ cannot be realized. 

1.1 Differences between NuTeV and R~ 

NuTeV does not measure exactly R~ but rather combinations of experimentally observed 
neutral-current to charged-current cross-section ratios, R meas in neutrino and anti-neutrino beams. 
Note that R~ itself can actually be written as 

R - ^(yv^ vo-^oyv^zvo 14 

a^N -> fx-X) - a^N -> p+X) ' 



aju^N^ii-X) ajV^N^ii+X) ajv^N-^^- X) 
a(v^N^ii-X) 

R v - rR v 



1 — r 



(1.5) 



(1.6) 



where R is the ratio of neutral-current to charged-current total cross-sections and r is the ra- 
tio of charged-current anti-neutrino to neutrino total cross- sections. Therefore, R~ itself is a 
combination of R from neutrino and anti-neutrino beams. 

NuTeV does not measure the total cross-section ratios of neutral and charged current inter- 
actions. The ratios actually measured by NuTeV, 

R^p = 0.3916 ± 0.0007 and if xp = 0.4050 ± 0.0016, (1.7) 



include non-muon neutrino backgrounds, the effects of experimental cuts, cross-talk between 
candidates in the numerator and denominator and final state effects. 

The electron neutrino background can be reliably subtracted as discussed below in Sec- 
tion and the effects of experimental cuts and cross-talk, although large, can be controlled 
experimentally and introduce no large corrections to the R~ dependence on sin 2 9w as will be 
shown in Section |4~21 

There are two final state effects which make modest corrections to the NuTeV observables. 
The first are small differences in acceptance in neutral and charged-current events due to the 
effect of the final state muon. The second is the effect of the final state in semi-leptonic charm 
decay. In such decays, some of the final state charmed hadron's energy is carried away in 
final state neutrinos and some may appear as final state muons. Further complicating the decay, 
those muons may affect whether the event is measured as a charged-current candidate or neutral- 
current candidate. These effects are small in the final analysis, but they do make some correction 
to the NuTeV observables not present in R~ . Both of these effects are fully modeled and 
corrected for in the NuTeV analysis. 

A more significant difference between NuTeV's analysis and R~ is the way that sin 2 6w 
is extracted from the two -R exp - The analysis uses a fit to two parameters, sin 2 ^ and the 
effective charm mass, m c , in charged-current charm production, which is the largest theoretical 
uncertainty in the analysis. This m c is externally constrained from the NuTeV measurements 
of charged-current charm production in the two muon final state, and so a 1C fit determines 
the result of Eqn. (I1.3I) . NuTeV has also performed a fit without this external charm mass 
constraint J24), and has found 

sin 2 #^ n " shell) = 0.22738 ± 0.00164(stat.) ± 0.00076(syst.), (1.8) 

which is consistent with the 1C fit result. 

It is instructive to decompose these (linear) fits into combinations of R^ xp and R^ xp in the 
spirit of Eqn. dl.61) . Writing 

b dsm 2 0w 

where the latter equality holds by choice of b, we find that for the 1C fit a = 0.2492, b = 0.6170, 
and for the OC fit a = 0.4526, b = 0.6116. The similarity between the b values for the two fits is 
a statement that the sensitivity to sin 2 9 W almost entirely resides in R v cxp and not -Rgxp' an ^ the 
near equivalence of the 1C and OC results indicates that R" xp is in agreement with expectations 
as is illustrated in Figure [U The discrepancy with the Standard Model of these results comes 
then from the fact that R" is not as expected. Note that these separate measurements, either 
cast as R v exp and R v exp or as the 1C and OC fits, will constrain any attempt to explain the NuTeV 
sin 2 9 W that seeks to make large changes in R u and R v separately to effect a large change in 
Br. 

2 Experimental Issues 

One of the largest experimental corrections to the NuTeV analysis is the subtraction of 
electron neutrino charged-current events from the neutral current candidate sample. Approxi- 
mately 5% of the neutral current candidates are electron neutrino charged currents [5 1, and so a 




Figure 1 : The measurements of R" xp and R" xp , shown as an error ellipse. The uncertainties 
in the error ellipse include theoretical uncertainties in relating i?g xp and R" xp to fundamental 
electroweak parameters, which result in the correlation between the two measurements. Note 
that the R" xp is in agreement with the Standard Model expectation, shown as the point, whereas 
the R" measurement is not. 



20% overestimate of this rate would be sufficient to explain the difference between the NuTeV 
sin 2 W and the Standard Model expectation. 

The dominant source of electron neutrinos in the NuTeV beams are decays. NuTeV 
determines the electron neutrino background by two methods: one an indirect determination 
using a beam Monte Carlo tuned to the neutrinos observed from decays, and the other a 
direct, but less precise, measurement of the electron neutrinos 0. An important check is that 
the direct and Monte Carlo methods agree, and in fact the ratio of measured to Monte Carlo 
predicted events is 1.05 ± 0.03 in the neutrino beam and 1.01 ± 0.04 in the anti-neutrino beam. 

The largest uncertainty in the Monte Carlo method is the 1.4% fractional uncertainty in 
the K% branching ratio @ . An interesting recent development comes from the BNL-E865 
experiment which has recently measured a branching ratio for K% |10| that is 6% larger than 
the value used by NuTeV. If this result is correct, it is interesting to note that it would not disrupt 
the agreement between the direct and Monte Carlo measurements of the u e rate at NuTeV and 
that it would in fact increase the discrepancy of the NuTeV sin 2 9w with the prediction by 
slightly less than one standard deviation ll24ll . 

The NuTeV observables, R" xp and -Rgxp' nave been studied as a function of position within 
the detector (Figures and 0), the containment length used to separate neutral and charged 
current candidates (Figure |H) and the visible energy in the detector (Figure Within the 
statistical and systematic uncertainties of these comparisons, no unexpected deviations as a 
function of event variables are observed. 

3 Electroweak Radiative Corrections 

NuTeV's analysis includes complete one-loop electroweak radiative corrections H2I3I . These 
corrections can be separated into those that can be absorbed into effective weak neutral current 
v-q couplings and those that cannot. The former cause a —0.00159 shift in the NuTeV reported 
sin 2 9 W . The latter corrections, dominated by a W-'y box diagram and by acceptance effects 
from bremsstrahlung in the final-state muon in the charged current shift sin 2 6w by an enormous 
-0.00795. 

One concern about the NuTeV result is that it relies on this single calculation of the radiative 
corrections. This calculation has been successfully checked against other partial calculations 
13 in the limits where both are expected to agree; however, another complete calculation as a 
cross-check might be advisable. 

4 QCD corrections 

A sizable part of the quoted theoretical error of sin 2 Q w as determined by NuTeV is due to 
QCD effects. Therefore, it is mandatory to investigate all potential sources of QCD uncertainties 
to the best possible accuracy. Sources of QCD effects include parton distributions and their 
uncertainties, perturbative QCD corrections at higher orders, isospin breaking effects or nuclear 
effects. In the following, we will discuss these in turn. 
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: The R cxp in the neutrino and anti-neutrino beams as a function of the depth of the 
interaction within the detector along the beam direction. 
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Figure 3: The NuTeV R exp binned in square annuli of transverse position from the center to the 
outer part of the detector. The first four bins are used in this analysis. 
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Figure 4: The effect of varying the NuTeV neutral-current/charged-current separation cut on 
R eX p in the neutrino (left) and anti-neutrino (right) beams. The error bars represent the sta- 
tistical uncertainty on the charge. The yellow band is the statistical uncertainty of the ratio as 
measured from the whole sample. 



4.1 Perturbative QCD corrections 

The QCD corrections to deep-inelastic scattering at higher orders are long known. For 
instance, in a complete next-to-leading order (NLO) treatment, and assuming massless quarks, 
the Paschos-Wolfenstein ratio Eqn. fll.ll ) becomes 1 



R- 



9l~9r + 



u 



dr 



(u- + d-) 2 



u + d~ 



3(#L - 9ru) + (did ~ 9m) + 7; — (9l ~ 9r) 

y 7T 



(4.1) 



Eqn. (14.11) uses the one-loop coefficient functions lfT3lfT4ll and holds, if no assumption on the 
parton content of the target is made. Here q~ = q — q is the second Mellin moment of the 
corresponding quark and anti-quark distribution, 



dxx (q(x) — q(x)) . 



(4.2) 



For consistency, NLO parton distributions are required, of course, and the nucleon is assumed 
to contain four light flavors. We will later comment on the treatment of charm as a massive 
quark. 

The result in Eqn. ( 14. II) has been expanded in powers of the dominant isoscalar combination 
of parton distributions, u~ + d~ . It shows the well-known fact, that the Paschos-Wolfenstein 
relation receives corrections, if the target has an isotriplet component, u ^ d, or sea quark con- 
tributions have a C-odd component, s~ ^ or c~ ^ 0. In particular, the QCD corrections only 

'In ref. [Ill the correct value [12| for SC 1 - 46C 3 is 32/9, with SC 1 = C 1 - C 2 and 8C 3 = C 2 - C 3 , cf. 
Eqn. O- 
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Figure 5: The R ex p as a function of visible energy in the neutrino and anti-neutrino beams. The 
ratio of data to Monte Carlo is shown below, with a green band to represent the systematic 
uncertainties in the comparison. 



affect these isotriplet u~ — d~ and C-odd terms, s~ ,c~. It is worth noting, that the coefficient 
at order a s in Eqn. (14.11) is the same in the M S-scheme and the DIS-scheme. This is due to the 
relevant combination of coefficient functions, 

- \c La + C 2 , q - C 3 , q (4.3) 

being invariant under changes between these two schemes. Of course, the NLO parton distri- 
butions entering in Eqn. (14.11) differ slightly in the two schemes, most notably in the region 
of larger x. Also, it is worth pointing out, that in the combination of coefficient functions in 
Eqn. (14.31) all dependence on the factorization scale /if (i.e. all logarithms \n(Q 2 / fi 2 )) cancels. 

Numerically, the correction factor u is —0.0232, evaluated using the NuTeV LO 

PDFs [5 1 at the experimental Q 2 which is extracted from fits to CCFR cross-sections lfT51 . In- 
clusion of the higher order terms in u -+ d - changes this to —0.0213. This evaluation assumes the 
strange and anti-strange seas carry equal momentum and that PDFs are isospin symmetric (see 
discussions in Sections l4~3l and l4~5l) . Therefore, the NLO QCD correction to R~ is —0.00035. 

It is possible, to extend Eqn. (14.11) to next-to-next-to-leading order (NNLO) in QCD, with 
the known two-loop coefficient functions lll61ll7[fT8lfT9ll . Putting the renormalization and the 
factorization scales /if = /i r = Q, we find 



d- 2 2 i ~~ r~ ° I o/ 2 2 \ , / 2 2 

R = 9l~9r + _ ; ,_ 3{g Lu - g Ru ) + {g Ld - g Rd 




1 , 16 , 83 1 a 2 

— G (q rif 

27 45 162 1 



Z§ \(9l ~ 9l) 



TV 2 



+ ° [ijFTFr j + ° y •) ' (4 ' 4) 



which holds in the MS-scheme. Here, rif denotes the number of light (massless) quark flavors, 
and C2 = 1.644934068, Cs = 1.202056903. For arbitrary scales fi f ^ Q, there will be at 
most single logarithmic dependence on Q 2 / /1 2 at order a 2 in Eqn. (I4.4I ). All double logarithms 
ln 2 (Q 2 //i 2 ) cancel due to Eqn. (|4.3b . In the DIS-scheme, the coefficient at order a 2 changes 
slightly. There, the expression in square brackets in Eqn. (14.41) becomes 



8704 5 . 16 . 83 



1215 9' 45" 162 



(4.5) 



Numerically, for typical values of a s , the NNLO corrections in both schemes, MS and DIS, are 
of the order of 30-40 % of the NLO contributions. This shows good perturbative stability of the 
QCD prediction with respect to higher orders. We have to note however, that a consistent NNLO 
analysis requires NNLO evolution of the corresponding parton distributions. The necessary 
three-loop anomalous dimensions are not completely known yet 12011211 . 

Thus, we can conclude the pQCD corrections to R~ are small and that uncertainties due to 
the perturbative expansion can be reasonably estimated. 

4.2 Experimental cuts 



Because of NuTeV's inability to measure NC events down to zero inelasticity, NuTeV cannot 
report ratios of total cross sections. One concern, then, is that R~ may not accurately reflect 



NLO QCD corrections to the NuTeV result. In particular, these cuts either remove or change 
the contributions near the kinematic endpoints in the inelasticity, y, which is a reason in general 
where one might expect enhanced sensitivity to radiative corrections. 

Backgrounds from sources other than muon neutrino interactions can be reliably subtracted 
in the NuTeV analysis; however, three major corrections to R~ remain. These are: (1) the 
minimum v cut required in order to observe the hadronic recoil in the detector, common to both 
the charged and neutral current, (2) the difference in the observed visible energy in the charged 
and neutral current due to the presence of the final state muon, and (3) charged-current events 
at very low which fake neutral currents. To a significant extent, all of these effects can be 
modeled by considering the relation dl.ll) with cuts on the inelasticity y. 

To that end, let us define for any difference of cross-section the cuts in y as 



{a(u.N^X)-a(V.N^X)} 



|/max 



2/min 



2/max 



dy 



da(u^N -> X) da(T7^N -> X) 



dy 



dy 



(4.6) 



with (standard) definition of the inelasticity y being the fraction of the neutrino's energy lost in 
the nucleon's rest frame. 

A simulation of the effect of the NuTeV broadband flux translates the experimental min- 
imum cut in visible calorimeter energy of 20 GeV into an effective minimum y-cut, y m i n , of 
0.24. Simulations of the energy deposited by the final state muon show a reduction of this ef- 
fective minimum y for charged-current events by 5 y of 0.03. Finally, a detailed detector and 
flux simulation Q can be used to measure the effective y at which the final state muons in 
charged-current events are so soft that the events are mistaken for neutral-currents. Denoting 
this y at which cross-talk occurs as 1 — y x , yx is numerically 0.043. 

This motivates the definition of a simple cross-section model for R~ as 



R 



model 



(ymin, Sy, yx) 



{a^N -> VO - a(v^N -> V^X)} ' + {a(^N -> n~X) - a{V»N -> /i+X)} 



(4.7) 
i 

l-ifx 



l-J/x 



which accounts for the kinematic cuts discussed. 

We can now work out the structure of -R mode i including higher order QCD corrections. The 
result for -R~ odel up to NLO can be written as follows (an extension to NNLO is straight for- 
ward): 

u~ — d~ + c~ — s~ 



^modci(2/min, 6y, yx) = /o(2/mi n , Sy, ^x) (flf| - g 2 R ) + hiVmin, Sy, yx) 

x (ju(ymin, Sy,yx)(gl U J till I 

+fd{ymin, yx)(gld ~ 9m) + /2(2/mm, Sy, y x ) 

+ — /a.(ymin, Sy, Vx){g L - g R ) + — h(ymm, Sy, y X ) 
7T 7T 



u + d 



1 



5 +Oiai,. 



\(u + d 

where the functions / , . . . , / 3> /„, f d and f as are given as 



(4.8) 



fo(vmin, Sy,yx) 
h(ymin,Sy,yx) 
fu(ymm,Sy,yx) 
fd(ymin,Sy,yx) 

f2(ymm,5y,y x ) 
fa s (ymin,Sy,yx) 



(jMm - l)(yLn - 2 ^/min ~ 2) 



-6 



P1P2 

(ymm ~ l)(Z/min ~ %min ~ 2) 
PlP? 

(ymin-l)(ymin-2y mi n-2) 
API 

x (?/min + y 2 x + Vx + yxSy + (1 + <5y) 2 - y min {2 + y x + 2Sy)) 

(j/min - 1 - $y)yx(-yx + y m in -1-5?/) 



-6- 



fsivmin, Sy,yx) 



g(Z/min l)(Z/min 2y m ; n 2) 



X ?/min - VnOnVX ~ ^minty ~ l% n 



-(13 + 5y)y x + 16 + y\ + 14cfy + <fy 2 

2 (y min - i - 5y)yx 
PiPi 



3 J l-i 2 

x (-3y x + 2y min y x - 7 + y min - Sy - 2y x Sy) , (4.9) 

where 

Pi = Z/min - 5y - 1 + y x 

P2 = yx-2- y min y x + y x - 2 2/mm + y 2 mhl + 25y + y x Sy + Sy 2 - 2y min Sy . (4. 10) 
In the limit y min , Sy, yx -> the functions / , . . . , / 3 , /«, fd and / a , simplify to 

foiVudu, Sy, yx) = 1, /i(l/mm, <fy, yx) = 0, 

fu(ymin, Sy, yx) = 3, fd{ymm, Sy, yx) = 1, 

fa s (y™in,$y,yx) 



9 



(4.11) 



in agreement with Eqn. d4.lt . Putting in the numerical values from the experimental y cuts, we 
find 

fo(ymm,5y,yx) = 1-053, fi(y m m, Sy,y x ) = 0.074, 
/ad/mm, <fy, yx) = 0.042, / 3 (y mm , <fy, yx) = 0.038, 
fu(ymm, Sy, y x ) = 2.700, fdiVwin, Sy, y x ) = 0.594, 

/a. (ymin, <fy, ?/x) = 0.683, (4.12) 



Numerically, then, the NLO QCD correction to -R^odei * s —0.00033. For comparison, the value 
for G?-R~ odel /cfem 2 9 W is —1.0075. In summary, in this model the experimental cuts in y do not 
result in large increases in the 0(a s ) QCD corrections. 

4.3 Parton Distribution Functions 

Let us next discuss the parton distributions. The very fact, that corrections to R~ in the QCD 
improved parton model are proportional to isotriplet or C-odd components of the nucleon target 
has led to questions about our knowledge of parton distributions for the various quark flavors. 
There are two major issues: the isotriplet component introduced by the small neutron excess of 
the NuTeV target, and possible momentum asymmetries between the strange and anti-strange 
seas. 

The NuTeV analysis corrects for the significant asymmetry of d and u quarks that arises 
because the NuTeV target, which is primarily composed of iron, has a SN = (A — 2Z)/A = 
+0.0574 ± 0.0002 fractional excess of neutrons over protons. As can be observed from the 
above uncertainty, the neutron excess is very well known due to a detailed material survey 
of the NuTeV target, including a chemical assay of the NuTeV steel Il22ll23ll24ll . The largest 
uncertainty in this isotriplet correction, in fact, comes from the difference between u and d 
PDFs. From the uncertainty on the NMC F^jF^ EH, NuTeV estimates this uncertainty to be 

0. 0003 in their extraction of sin 2 9w- (Note that this correction assumes isospin symmetry, i.e., 
(-) / \ ( } . . ( ) . ^ (-) , 

u p(x) = d n{x), d p{x) = u n (x), violations of which are discussed below.) 

Kulagin [ 26 ] has recently noted that NuTeV's measurement of the neutron excess correction 
in sin 2 9 W , 0.0080, does not agree with the correction of the effect in R~, and suggests that the 
difference could be taken as a (substantial) correction to the NuTeV measurement. However, the 
differences can be partly understood in terms of the experimental cuts discussed in Section |4~2~1 
which reduce the neutron excess correction by 20%, and partly understood as the differences 
between the NuTeV 1C fit and R . Therefore, the suggestion by Kulagin that this difference is 
a correction to NuTeV's sin 2 9w is incorrect. 

Global analyses of unpolarized parton distributions usually assume no strange asymmetry, 

1. e. imposing as a constraint s(x) = s(x). Moreover, a fit ansatz like ll27ll 

/x / \ u(x) + d(x) ,,,..„* 
s(x) = s(x) — K K ' 2 K J , (4.13) 

ties the strangeness distribution to the relatively well known u(x) and d(x) distributions, thereby 
underestimating the true uncertainty [28|. There are several recent dedicated analyses of inclu- 
sive lepton-nucleon DIS (including neutrino data) I29II30II31L which have used various methods 
to account for the uncertainties and correlated errors in parton distributions [ 30 1, ll32l . 

One of these analyses of the CDHS neutrino data and charged-lepton data reported 
some improvement in their fits if they allow for an asymmetry in the strange sea at high x. 
However, this large asymmetry at high x is directly excluded by the non-observation of high 
x dimuon events at NuTeV and CCFR ll33ll . An update of this analysis using CCFR and 
CDHS charged-current neutrino data no longer finds the large significant asymmetry between 
the strange and anti- strange quark distributions lf34l . 



As illustrated above, the NuTeV dimuon data 031l can help in understanding the quark 
flavor content of the nucleon as it provides separated measurements of s(x) and s(x) initial 
state contributions to charm production. 

NuTeV has measured possible differences between s(x) and s(x) from the CCFR and 
NuTeV data on un^n —> n + n~X within the NuTeV enhanced LO cross-section model used 
in the sin 2 9 W analysis. Denote the momentum integrals of the strange and anti-strange seas as 
S, S, i.e., S = J xs(x)dx and S = J xs(x)dx. Within this model, the dimuon data implies a 
negative asymmetry |22J, 

S- S = -0.0027 ±0.0013, (4.14) 

or an asymmetry of 11 ± 6% of (S + S). Therefore, dropping the assumption of strange- 
antistrange symmetry results in an increase in the NuTeV value of sin 2 9w, 

Asin 2 9 W = +0.0020 ± 0.0009. (4.15) 

The initial NuTeV measurement, which assumes s(x) = s(x), becomes 

sin 2 ^ n " shcll) = 0.2297 ± 0.0019. (4. 16) 

A preliminary analysis of the strange and anti-strange asymmetry in an NLO cross-section 
model also finds the momentum carried by the seas to be consistent within uncertainties. It is 
worth noting that these fits have been criticized in the literature [ 1 1 1 because of their assumed 
functional form in x, and indeed, NuTeV is planning to reanalyze this data with functional 
forms more suggestive of strange and anti-strange asymmetries suggested by non-perturbative 
calculations IE5TOI571I5SI. 

The impact of uncertainties in the strange-antistrange and in the up-down quark asymmetries 
of the target on the NuTeV measurement of sin 2 9 W has been quantified in terms of an error 
functional E3, F[sin 2 9 w ,5;x] such that 

Asm 2 9 w =[ Flsin 2 9 w ,S;x}S(x) dx, (4.17) 
Jo 

for any symmetry violation S(x) in PDFs. All of the details of the NuTeV analysis are included 
in the numerical evaluation of the functionals shown in Figure For this analysis, it can be 
seen that the level of isospin violation required to shift the sin 2 0w measured by NuTeV to 
its standard model expectation would be, e.g., / x(d p (x) — u n (x))dx ~ 0.01 (about 5% of 
/ x(d p (x) + u n (x))dx), and that the level of asymmetry in the strange sea required would be 
S - S ~ +0.007 (about 30% of S + S). 

4.4 Charged current charm production 

Thus far, the discussion of QCD corrections to R~ as in Eqns. (14.11 1. (14.41) has assumed 
quarks to be massless. However, an important contribution to the charged current reactions is 
the deep-inelastic charm production in us — ► fi~cX. Clearly, the mass m c of the charm quark 
should to be taken into account. 

For the center-of-mass energies accessible to NuTeV, a scheme of three light flavors in the 
nucleon is expected to give an appropriate approximation. Thus, charm quarks are entirely 
generated perturbatively. The leading order QCD improved parton model then requires the 




Figure 6: The junctionals describing the shift in the NuTeV sin 2 6w caused by not correcting 
the NuTeV analysis for isospin violating u and d valence and sea distributions or for (s(x)) ^ 
(s(x)). The shift in sin 2 9w is determined by convolving the asymmetric momentum distribution 
with the plotted functional. 



x-dependence to change slightly. The relevant strange distribution becomes s(£, //f), where 
£ = x(l + m 2 /Q 2 ) as the scaling variable correctly accounts for the single-charm threshold 
condition. This amounts to the so-called "slow rescaling" [39|. This, for instance, has been 
implemented in the NuTeV analysis of the dimuon production cross-section ll33l . 

Also the higher order QCD corrections to us — > fj,~cX are known since long. At or- 
der a s , the complete corrections to deep-inelastic charged current scattering have been calcu- 
lated H391I40I . From these results, it is straight forward to obtain the complete m c -dependence 
of R~ at NLO. However, for the purpose of this letter we restrict ourselves to a discussion of 
certain qualitative features. First of all, for an observable like R~, consisting of a particular 
combination of total cross-sections, some dependence on charm mass cancels. Most promi- 
nently, there will be no large logarithms ln(Q 2 /m 2 ) in the NLO corrections to R~ due to the 
combination of coefficient functions in Eqn. (14.3b . In comparison with Eqn. (14.11 ). the treatment 
of charm as a massive quark will at most introduce additional terms of order m 2 /Q 2 in the 
coefficients proportional to a s . 

At order a 2 , only those corrections for deep-inelastic charged current scattering are known, 
which are logarithmically enhanced by 0(a 2 ln n (Q 2 /m 2 )) contributions BT1 . For R~ this 
implies, that there can be at most single logarithms \\i(Q 2 /m 2 ) at NNLO due to the absence of 
logarithmic enhanced terms at NLO. 

Thus, we conclude that for R~, any dependence on the charm mass should be weak. The 
same conclusion also holds, if we model the effect of experimental cuts like in -R^odei °f 
Eqn. (14.71) . There again, the m c -dependence enters only through additional terms of order 
m 2 JQ 2 . 

However, as already mentioned in Section fTTTl the charm mass dependence does play a role 
in choosing how the experimental R u and R v are combined to extract sin 2 8 W . Although the 
result is robust over changes in this prescription, this dependence is still worth noting and may 
bear further investigation. Also, the extraction of the strange sea from the CCFR/NuTeV dimuon 
data, which is used as input to the sin 2 6\y measurement has significant NLO corrections. Again, 
the (quark-antiquark symmetric) strange sea doesn't enter directly into R~ but does enter into 
the experimental combinations used in the NuTeV fits. 

For a full quantitative analysis of the various experimental cuts in phase space a detailed 
Monte Carlo study has to be performed. This has to account also for fragmentation, which 
has been entirely neglected in this discussion. In particular, the fragmentation has a non-trivial 
effect on the visible final state energy due to the presence of the hard neutrino from the charm 
decay. A Monte Carlo program, which calculates the fully differential cross-sections has been 
provided in 11421 . 

4.5 Isospin breaking 

Having concluded the investigation of perturbative QCD corrections to R~ let us now turn 
to investigate potential nuclear effects and isospin violation. 

As argued above, knowledge of the neutron excess allows for a reasonably accurate cor- 
rection for the isovector part of the cross-section; however, this correction is only valid with 

(-) , x H . . (-) . . (-) . . 
the assumption of isospin symmetry, i.e., u p (x) = d n \x), d P {x) = u n (x). This as- 
sumption, if significantly incorrect, could produce a sizable effect in the NuTeV extraction 

of sin 2 e w MMMMMMQ- 



Let us briefly discuss the main proposed non-perturbative models to generate isospin viola- 
tion in the nucleon I1431l44ll431l . The earliest estimate in the literature, a bag model calculation 
by Sather [43 1, predicts large valence asymmetries of opposite sign in u p — d n and d p — u n at all 
x, but neglects a number of effects. Most notably, the effective mass of the remnant diquark is 
assumed to have a 5-function distribution. Recently, Londgergan and Thomas, revisiting their 
earlier calculation [44] with the fixed diquark mass in the bag model but including a number of 
effects neglected by Sather, including nucleon size and mass, have argued that their calculation 
observes the same effect as Sather [46 ] and that this effect is largely independent of PDFs PP. 
at least when the x dependence of the NuTeV acceptance ll22ll is neglected. However, when 
this same calculation is done with a smeared distribution of diquark masses H4l . the dominant 
isospin violating effect of the minority quark distribution, d p — u n , is reduced at high x and the 
negative asymmetry at low x is found to carry more momentum. Thus including the effect of di- 
quark smearing, the high x and low x contributions largely cancel, leaving a negligible (0.0001) 
shift in the NuTeV sin 2 9w E3- The effect is also evaluated in the meson cloud model [45|, and 
there the asymmetries are much smaller at all x, again resulting only in a small shift in sin 2 9w- 

Models aside, the NuTeV data itself cannot provide a significant independent constraint on 
this form of isospin violation. However, if parton distributions extracted from neutrino data (on 
heavy targets) are used to separate sea and valence quark distributions which affect observables 
at hadron colliders ll48ll . effects of isospin breaking could be seen. Therefore, global analyses of 
parton distributions including the possibility of isospin violation may be able to constrain this 
possibility further experimentally. 

4.6 Nuclear effects 

Any nuclear effect which can be absorbed into process-independent PDFs will not affect 
the NuTeV result. However, several authors have recently explored the possibility that neutrino 
neutral and charged current reactions may see different nuclear effects and therefore influence 
the NuTeV result. 

Thomas and Miller [ 49 ] have offered a Vector Meson Dominance (VMD) model of low x 
shadowing in which such an effect might arise. The NuTeV analysis, which uses v and V data 
at < Q 2 > of 25 and 16 GeV 2 , respectively, is far away from the VMD regime, and the effect 
of this VMD model is significantly smaller than stated in this analysis. The most serious flaw 
in the hypothesis that this accounts for the NuTeV result, however, is that it is not internally 
consistent with the NuTeV data. Shadowing, a low x phenomenon, largely affects the sea quark 
distributions which are common between v and V cross-sections, and therefore cancel in Rr . 
However, the effects in R u and R v individually are much larger than in R~ and this model 
increases the prediction for NuTeV's R v and R v by 0.6% and 1.2%, respectively. NuTeV's R v 
and R v are both below predictions and the significant discrepancy is in the v mode, not the V 
"control" sample, both in serious contradiction with the prediction of the VMD model ll50l . 

Kulagin [26| has recently investigated possibilities for process-dependent nuclear effects 
disrupting NuTeV due to Fermi motion, nuclear binding corrections and shadowing, and found 
the effects to be small. Schmidt and collaborators [5 1 1 have suggested that there may be little or 
no EMC effect in the neutrino charged current but the expected EMC effect suppression at high 
x in the neutral current. If true, this could have the right behavior and perhaps magnitude to 
explain the NuTeV data because of the effect at high x. Unfortunately, this mechanism would 



1.3 
1.2 
1.1 
1.0 
1.3 
1.2 
1.1 
1.0 

1.20 
1.15 
1.10 
1.05 
1.00 
0.95 
1.10 

1.05 

1.00 

0.95 

1.10 

1.05 

1.00 

0.95 

1.10 

1.05 

1.00 

0.95 





0.5 1.0 5.0 10.0 50.QL00.0 

2 



0.5 1.0 5.0 10.0 50.QL0O.O 

n2 



TR-VFS (MRST99) |>=Q] 

- - - - ACOT-VFS (CTEQ4HQ) [^=AC|OT] 

- ■ - FFS (GRV94) [/i=2mc] 
NMC o BCDMS o SLAC 



Figure 7: The ratio of F% with a model-independent prediction from F| on iron. Primarily 
because of quark-mass effects these ratios are not expected to be one, and different theoretical 
calculations for this ratio are shown. Good agreement of the ratio with expectations limits 
anomalous nuclear effects in the v charged-current interaction. 



cause large differences between F% and F| on heavy targets at high x which are excluded by 
the CCFR charged current cross-section measurements lfT51l shown in Figure |7| 

Kumano [52| has fit experimental DIS and Drell-Yan data on nuclear targets to investigate 
the possibility that nuclear effects are flavor dependent, rather than process dependent. Such 
an effect could impact NuTeV because the constraints on d/u of the nucleon come from light 
targets. Kumano performs these fits in the context of "nuclear PDFs", and found very small 
flavor-dependent effects, except at very high x and low Q 2 , a region removed by the visible 
energy requirement (^calorimeter > 20 GeV) of the NuTeV analysis. The effect on the NuTeV 
analysis which assumes flavor-independent nuclear effects is therefore negligible. 

5 Conclusions 

Motivated by the NuTeV measurement of sin 2 9w, which deviates from the standard model 
prediction by approximately 3 a, we have studied potential sources of theoretical uncertainties. 
In particular, we have investigated QCD effects at higher orders in perturbation theory including 
the dependence on the charm mass. We have discussed the impact of parton distributions and 
potential nuclear effects. To assess the effect of experimental cuts on the measured R~, we have 
developed a simple model -R~ odel which accounts for y-cuts up to NLO perturbative QCD. 

Based on these investigations, we conclude that higher QCD corrections are under con- 
trol and small. The uncertainties on parton distributions have been discussed and will also be 
addressed in future global analyses. An asymmetric strange sea seems unlikely to be an ex- 
planation of the present discrepancy. Large isospin violation in parton distribution functions 
is a possible explanation, but the violation would have to be larger than naive estimates would 
suggest. 
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Here we give some formulae relevant for the calculation of QCD corrections to R . Neu- 
tral/charged current cross-sections are defined as [9| 



Appendix 
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dxdy 



(5.1) 



with 



Y + 



(l + (l-2/) 2 )/2 



Y- = {l-{l-yf)/2. 



(5.2) 



F? c (x) = x 



The k NC / cc in Eqn. (15.11) are expressed in terms of the structure functions F 2 , F 3 , F L 

fc NC, = Y + Ff c + Y-xFf c -y 2 /2F* c , 

k NCu = F + F NC_y- a . jF NC_ y 2 /2F NC 

k ccu = Y + p cc + y- xF cc _ y i/2F£ c , 

k ccu = Y + F cc _ Y ~ xF cc _ y 2/2F£ c . (5.3) 

For the Paschos-Wolfenstein R~, the discussion can be restricted to the non-singlet structure 
functions. In the QCD improved parton model (assuming massless quarks), these are given as 
convolutions of parton distributions and the coefficient functions of the hard scattering process 
(see for instance H 14111 8H or Handbook of perturbative QCD J53J). 

' {(u 2 L + u 2 R )(u(z)+u(z) + c(z) + c(z)) 

+ (4 + d 2 R )(d(z) + d(z) + s(z) + s(z))]C iiq (x/z) , 

(x) = x [ — Uu\ - u 2 R ){u{z) - u(z) + c(z) - c(z)) 

J x Z K 

+ (4 - d 2 R ){d{z) - d(z) + s{z) - s(z))}C 3a (x/z) , 

rl 

F l cc (x) = x — {u(z) + d(z) + s(z) + c(z)}C i>q (x/z) , 

J x Z 

xF 3 cc (x) = x — {-u(z) + d(z) + s(z) -c(z)}C 3a (x/z) , 

J x Z 

F l cc (x) = xj y [u{z) + d(z) + s(z) + c(z)]C in (x/z) , 
xFp(x) = x f 1 — lu(z) - d(z) - s(z) + c(z)\C 3 , q (x/z) , (5.4) 

J x Z K 

where % = 2, L and the non-singlet coefficient functions have an expansion in powers of a s , 

„2 



xF 3 



NC/ 



C 2 , q (x) = 5(1 - x) + -^c^(x) + — ^ >q (x) + 
C 3 , q (x) = 5(1 - x) + ^cgj(x) + j^cflix) + 



Cl«{x) = ^c«(x) + ^ci 2) q (x) + . . . (5.5) 

If restricted to leading order, these equations reproduce the naive quark parton model. 

In Eqns. (15.41) . (15.51) all scale dependence has been suppressed. The structure functions, be- 
ing observables, depend on the physical scale Q, whereas the PDFs and the coefficient functions 
depend on the renormalization scale fi T and the factorization scale /if, the usual choice being 
fji = fir = fif. Putting additionally /i = Q also cancels all scale dependent logarithms in the 
coefficient functions at Woops, i.e. all terms ln*(Q 2 / /j, 2 ) , i < I. 

In the Paschos-Wolfenstein relation R~ only total cross-sections enter. Here, the expressions 
in Eqn. (15.41) simplify considerably. After the x-integration to obtain total cross-sections, all 
convolutions become simple products of (Mellin)-moments. Thus, 



J dxFf c = {(u 2 L + u 2 R )(u + u + c + c) 



+ (d 2 L + d 2 R )(d + d + s + s)}C iA , (5.6) 

J dxxF^ c (x) = {(u 2 L -u R )(u-u + c-c) 

+(d 2 L -d 2 R )(d-d + s-s)} C 3 , q , (5.7) 

f 1 dxFp c (x) = {u + d + s + c}C iq , (5.8) 
Jo 

[ dxxF^ c (x) = {-u + d + s -c}C 3 q, (5.9) 
Jo 

where i = 2, L and the quantities on the right hand sides are the second Mellin moments of 
PDFs and coefficient functions. Beyond leading order, the factorization of structure functions 
into PDFs and coefficient functions is arbitrary and introduces scheme dependence. The most 
common schemes are MS, for which most higher order cross-sections have been calculated, 
and the DIS-scheme, which is physically motivated by demanding 

C 2 , q (a;) = 5(1 -x) (5.10) 

to all orders in perturbative QCD. Of course, the PDFs and the coefficient functions change 
accordingly and the implications for the Paschos-Wolfenstein relation have been discussed in 
Section 4. 
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